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Abstract 

In this paper we compare two classical one-factor diffusion mod- 
els which are used to model the term structure of interest rates. One 
of them is based on the Wiener-Bachelier process while the second one 
is based on the Ornstein-Uhlenbeck process. We show essential differ- 
ences between the prices of European call options on a zero-coupon 
bond in these models. 

PACS numbers: 02.50.-r, 03.67.-a 

Keywords: Stochastic processes, Merton model, Vasicek model, Wieiener- 
Bachelier process, Ornstein-Uhlenbeck process 



1 



Introduction 



The bond market forms a very important segment of financial markets. 
However, modelling of it is difficult because the bond prices depend on both 
the present situation of the market and the time to maturity. The prices of 
bonds are expressed in terms of various interest rates and yields, so under- 
standing bonds pricing is equivalent to understanding interest rate behaviour. 
The term structure of interest rates is defined as the dependence between in- 
terest rates and maturities. We propose two models: 

1) the Merton model based on the Wiener-Bachelier process, 

2) the Vasicek model based on the quantum Ornstein-Uhlenbeck pro- 

*T)\ 
cess|__: 

In the above models we find the formula for the zero-coupon bond price 
and we price the option on this bond. In the first section we propose gen- 
eral zero-coupon bond pricing model and we give some information about 
forward contracts with the help of which we price options on zero-coupon 
bonds. Then, we discuss the Merton model while in the third one we discuss 
the Vasicek model. Finally, we compare both models. 



We assume that the market is effective without an opportunity of arbi- 
trage and the face value of any bond equals one. Let Bf be the price of 
the zero-coupon bond at some date t with the time to maturity T > t. Let 
yj denotes the discount rate called zero-coupon rate or spot rate for date 
T in continuous capitalisation [S]- Then Bj is related to the continuously 
compounded zero-coupon rate by 



hence 



yl=-—-\nBi. (2) 



We call the function T — > Bj the discount function and the transformation 
T — ^ y'^ the yield curve. 



The short-term interest rate is given by the formula 

rt = lim yl'^'^ ■ 

r— »0+ 

*^)ln the quantum game theory Ornstein-Uhlenbeck process has the interpretation of 
non- unitary tactics resulting in a new strategy. These strategies are the Hilbert's spaces 
vectors. Quantum strategies create unique opportunities for making profits during intervals 
shorter than the characteristic thresholds for the Brown particle, see [TJ [2] . 
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We introduce few notions which are connected with the derivatives. First, 
we consider options on the zero-coupon bonds. Let us consider European 
options. We must clearly distinguish between the call option and the put 
option. The time of maturity of the option is denoted by T. Let Ct denotes 
the time t price of an European call option on the zero-coupon bond, which 
gives a payment of 1 at time S, where S > T, with the exercise price of 
K. Similarly, let ttj denotes the price of the corresponding put option. The 
prices of these options at maturity are equal to Ct = max(i?|i — K, 0) and 
ttt = max{K - S|, 00. The prices of European call and put option on the 
zero-coupon bonds fulfil the put-call parity relation [4] 

Ct + KBj = + fif . (3) 

As the result, we can restrict ourself to pricing of European call options. 

Another derivatives are forward contracts on the zero-coupon bonds. If 
the underlying variable is the price of the zero-coupon bond, the essential 
part by the pricing of the forward contracts plays the present prices of the 

T S 

zero-coupon bonds only. If there is no arbitrage, the time t value ' of the 
forward contract on the zero-coupon bond, with the delivery date T and the 
delivery price K, is given by the formula [4j 

y,^'^ = Bf - KBj. (4) 

5" > T is the maturity date of the underlying bond. For forwards contracted 
upon at time t, the delivery price K is set so that the value of the forward 

'T S T S 

at time t is zero (V^ ' =0). This value of K is called the forward price ' . 
Solving (jll) we obtain that 

= (5) 



1 The pricing in affine models 

We assume that for considered market model, there is a risk-neutral prob- 
ability measure (or equivalent martingale measure) P* and a one-dimensional 
Wiener process W* under measure P*. If the risk-neutral probability measure 
exists, then there is also the T-forward martingale measure P^ and then W'^ 
is the Wiener process under measure P^ . If there exists a spot martingale 
measure, then the market is arbitrage free. Part of the market characteris- 
tic fulfils the short rate under the risk-neutral probability measure (i.e. the 

*'^^ BJ always denotes the price of the zero-coupon bond at the time t with the time to 
maturity T . 
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spot martingale measure). The model in which this rate is represented by the 
process {rt)t>o, that fulfils the stochastic equation 

dn = fi{rt, t)dt + cr(ri, t)dW^ , (6) 

where ^ : M x [0,T] — ^ M and a : R x [0, T] — > W^, we call the affine 
model. Let us consider that Ht = Oii — a2 ■ n and at = y/ Pi + ■ rt, where 
ai, a2, Pi, P2 are constant. We know that the discount price of the zero- 
coupon bond is given by the formula [3] 



Bi =E 



p* 



(7) 



With the help of this assumption, we can calculate the formula for the dis- 
count price of the zero-coupon bond, see This price can be written as a 
function of time and the current short rate r. 

Theorem 1.1. 

In an affine model, in which the short rate r is described by the time t 
price of the zero-coupon with the time to maturity T and face value 1 is equal 
to 

= ^-a{T-t)-rKT-t) _ ^g) 

The functions a,b : [0,T] — > M fulfil the following system of ordinary differ- 
ential equations: 

^P2b{ty + a2b{t) + b'{t) -1 = 0, 

a'{t) - aib{t) + ]^(3ib{tf = 0, 

for t G (0, T) and a(0) = 6(0) = 0. 
Proof. 

The proof is based on the Feynman-Kac theorem, see flS.Sp in Appendix. 
From theorem (15.81) follows that Bj, given by the formula ([H]), fulfils the 
partial differential equation 

^ + («i - a2r)^ + -iPi+ P2r)^ - rBf = 0. 

The relevant derivatives are 
dB7 



't _ TjT 



dt 



Bf{a'{T-t) + rb'{T-t)), 



= -BibiT-t), 

f)2 tdT 

= Bjt,[T-tf. 
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After substituting these formulas into fl5.8p . we obtain 



a\T - t) + rb\T - t) - (ai - a2r)h{T - t) + h{T - tf{l3^ + /^sr)^ = 0. 
After gathering terms involving r, we obtain 

(a'(T - t) - ai6(T - t) + ^/5i6(T - + 

+ (i/?26(r - if + asK^ - + ^'(^ - ^) - 1)^ = 0. 

The last equation must be fulfilled for arbitrary r, so the expressions in the 
brackets must be identically 0. Hence, we find equations which the function 
a and h must fulfil. Because = 1, so a(0) + 6(0)r = and hence a(0) = 
6(0) = 0. We see that the function fulfils the Feynman-Kac equation. 

□ 

Let us observe that, if we have the function b from the above formula, we 



can calculate the function a from the formula*^-* 

a{t) = ai [ b{s)ds - [ b{sfds. (9) 
Jo 2 Jo 

If we want to determine the price of the bond, we have to find the solution 
of the differential equation given by the formula 

\l32b{tf + a2b{t) + b'{t) -1 = 0, 6(0) = 0. 

Now, let us observe that bif) > for t > 0. Above all, from the continuity 
of b, the conditions 6(0) = and b'{t) = 1 - \(32b{tf - a2b{t) follow that 
b'{t) > in the vicinity of 0. From this, b grows up in the vicinity of 0, that 
is b{t) > for small t > 0. As if b{t) < for certain t > 0, there is a point 
ti > from the Intermediate Value Theorem [7j, such that 6(ti) = 0. From 
the Rolle theorem [7], we obtain 6(0) = 6(ti) = 0, so there must exist to > 
such that 6'(to) = 0. If we substitute to to the differential equation fulfilled by 
6, we obtain the equality |/32?/o + '^^Uo — 1 = for ?/o = 6(to)- Let us consider 
constant function 6o(t) = yo. It fulfils the same differential equation, so 6. 
Furthermore, 6(to) = 6o(to); that means that both functions fulfil the same 
initial-value problem. From the uniqueness of solution we obtain that 6 = 6o 



*3)lf a'{t) = aib(t) - ^Pib{t)^, so a{t) = a(0) + ai b{s)ds - /J b{sfds. Let us take 
into consideration that a(0) = 0. 
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and 60(0) = f'(O), that is, yo = 0, which is impossible because does not fulfil 
the quadratic equation which is satisfied by yo. We obtain contradiction, so 
b{t) > for t > 0. It means that the function r — > Bj is strictly decreasing. 

In the further parts of our paper we will use the formula for the volatility 
of the process under T-forward martingale measure. It looks the same as 
for measure P*. From the Feynman-Kac theorem, it is given by the formula 



9BT 



and in our specific situation fit = y/ (ii + /?2r. Using the formula ([H]), we obtain 

= '^^^J^^' Pt = -h{T - m. (10) 

Let /i^ denotes the relative drift of the process B^ under measure . 



2 The Merton model 

The first dynamic, affine and continuous time model of the term structure 
of the interest rate was described by Merton |8j. The short rate follows a 
generalised Brownian motion under the spot martingale measure: 



drt = ^pdt + adWf 



t ) 



where a and if are constant. We begin by finding the price of the zero-coupon 
bond and then we calculate the price of the option. We see that the price of 
the bond in this model is given by the formula (see ([H])) 

where the function b fulfils the simple ordinary differential equation b'{t) = 1 
with 6(0) = 0. So, b{t) = t and the function a is given by the formula (see 

©) 

/"* 1 /"* 11 

ait) = tp sds (T^ / s'^ds = -iot"^ a^t^, 

Bj = e-|¥'(T-t)2 + ia2(T-t)3-(T-t)r^ (^^^^ 

In order to price the option, we use the theorem (15.71) . see Appendix. The 
price of an European call option on the zero-coupon bond with the expiration 
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date S, the exercise price K, and maturing at the time T is given by, see ffT^ 
in Appendix, 



Ct = Bj (e^^ [max(S| - K^^t]) (r). 



(12) 



From ([5]) we know that the forward price fulfils the formula: 



Bl 



hence, in particular F^'^ = B^. The forward price (F7'''^)f>o is a martingale 
under the T-forward martingale measure, what means that the drift rate 
equals under this probability measure. Simultaneously, this process is a 
quotient of the prices of two zero coupon bonds. In the diffusion equations for 
B^ and B^ , the relative volatilities of the bond are equal to erf = —a - {T — t) 
and erf = —a ■ {S — t) respectively, so that 



dB'l' = ijjBjdt - a{T - t)B'i'dWt 



t 1 



dB^ = fifBtdt - a{S - t)B^dWl . 

By an application of the Ito lemma, see Appendix (15.31) . for functions of 
multiple stochastic processes and because of that, the drift of the process 
pT,s gquaig g under the T-forward martingale measure, we obtain 



dR 



T,S 



dB 



T t 



dBf 



dw; 



B 



(tTbT + 



dw: 



= (af - aj) Fj-'dWf = ~a{S - T)Fr'dW^. 

This means that, F^'^ follows a geometric Brownian motion, see Appendix 
(15.41) . with a'^'^{t) = —a ■ (S — T). Hence, the random variable X := In 
is normally distributed and 



X 



a'^'^{uYdu + 



a 



[u)dWff, 



(13) 



what result from (15.51) . see Appendix. Moreover, 



E{X) 



a^'^{uydu = -^a\S-T)\T-t), 



v{t, T, S) := ^/Var{X) = J a'^^^iufdu = a(^ - T)VT^t. 



Let us consider the function (7 : M x M_,_ 3 (x, y) — > ma.x{e^y — K, 0) and a 
random variable Z = . From the formula (fT^ and with the help of that 



= F^'^ , we obtain 



Ct =BT E 



max(F^'^ - K,0)\J^t 



From the equation ( fT3l) we have that X is independent of J-'t and Fp^ is 
measurable, because the process F^''^ is a martingale. Hence, we can use the 
lemma g^H see Appeudix. By the using the lemma (15.61) . see Appendix, 
we can calculate the expected value [?/ > !] 



E{g{X, y)) = E{max{Z ■y~K,0)) = yE{ max 



^E{X)+\v(t,T,Sf ^ 



fE(X)-\n^ 



-N 



vit,T,S) 



.eOiv(-l.(t,T,^) + -^ln| + .(t,T, 



S) 



In^ 



yN 



^ ^n^ + lv{t,T,S)]-KN 



v{t,T,S) K 2 



1 V I 

■In^ --v{t,T, S) 



v{t,T,S) K 2 



T S 

If we substitute to the above equation y = F^ ' , we obtain 



E 



g{X,Fl^')\J^, 



KN 



In 



v{t,T,S) K 2 



and finally 



c, = bT 



F^'^^Nidi) - KN{d2) = B^N{di) - KB'l'N{d2 



T ] 



**'>Let us observe that g{X,Ft ^) = inax{Fr^'^ - K,0), so that g{X,F'^^^) is bounded, 
because < F'^^'^ < 1. 
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where 



v{t, T, S) = a{S - T)VT^. 

3 The Vasicek Model 

In the paper [H] Vasicek proposed a mean-reverting version of the Ornstein- 
Uhlenbeck process for the short term rate. Specifically, under the risk-neutral 
measure P*, rt is given by 

drt = k{9 - rt)dt + adW;, 

where k, 6 and a are positive constants. The parameter 6 denotes long-term 
level of the short rate, because the rate vt pull towards a long-term level of 
9, K determines the speed of adjustment, and a is the average deviation of 
the rate of return. 

Similarly to the Merton model, at the first we calculate the formula for 
the zero-coupon bond price. From the theorem 11.11 we know that 

^ ^-a(T-t)-rb{T-t) ^ 

but in the described model we have 

Kh{t) + b\t)-l = Q, 6(0) = 0. 

Hence, 

h{t) = 1(1-6-''*). 
K 

From ([9]) we obtain that 



a{t) = kO f b{s)ds - [ b{s)'^ds 
Jo 2 Jq 

2k' Jo 



t 2 ft 

9 I {1- e~''')ds j (1 - 2e-"* + e-^''')ds 
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If we substitute a{t) and b{t) to the formula ([H]), we obtain analogous formula 
for the bond price to the Merton model, see (fTTj) . 



We calculate the formula for the call option price similar to the Merton 
model. Let us use the basic formula 

Ct = BJE^^ [max(S| - ^,0)1^;] . 

With the help of the Ito lemma f l5.3p and because of that the drift of pro- 
cess F'^'^ equals under the T-forward martingale measure, we obtain the 
diffusion equation for the forward price: 

dF^'' = ((af - aj) Fj^'dWj = -a [h{S - t) - b{T - t)] F]^''dW^. 

Therefore, the process of forward prices is a geometric Brownian motion with 



T,S 



(t) = —a \b{S — t) — b{T — t)]. The random variable X := In has the 



Gaussian distribution. Similarly, we obtain the formulas: 




u)dWl 



Our further reasoning is the same like in the Merton model. Let the function 



g : Kx 
We obtain the formula 

Ct 



3 {x, y) I — > max(e^y — K, 0) and the random variable Z = e 



X 



T S 

With the help of independence X of J-'t and J-'t-measurability of F^ ' , we can 
use the lemma (15.21) . Then, if we use the lemma fl5.6l) . we calculate that the 
expected value E{g{X,y)) equals 



-KN 



-^^in|-i.(t,r,5)^ 
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T S 

If we substitute to the above equation y = ' we obtain 



1 



pT,S 



and finally 
where 



KN 



In 



v{t,T,S) K 2 
Ct^BfN{d,)-KBjN{d2), 



7^v{t,T,S) 



di 
d2 

v{t,T,S) 



1 



In 



Si 



1 [ \ 1 



(1 - e-''[^-*l)Vl-e-2«[^-*l. 



4 Final remarks 

The Fig. 1 shows the prices of European call options on the zero-coupon 

bonds with the expiration date 5 years, the face value 1 and the real val- 
ues of the market parameters based on the Merton model and the Vasicek 
model. In the above figures, we consider the European call option on the 




Figure 1: The prices of European call options on the zero-coupon bonds based 
on the: 1.1- Merton model and 1.2-Vasicek model. We assume that: k = 0.4, 
a = 0.03, r = 0. 
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zero-coupon bond with the exercise price K = 0.80, time to maturity 5 = 5 
years and the face value i?f = 1. As it follows from the general affine model 
pricing, see (E]), we can compare this two models under the assumption that 
(f = k9. The results are qualitatively different. Let us notice that, the price of 
the option described by Vasicek is higher than the price described by Merton 
and faster drawn towards to 0.2, where 0.2 is the expected price of the option 
with fixed expiration date T = 5, because the market is effective without an 
opportunity of arbitrage. Let us observe that, both models give equal prices 
for the long interest rate 6* = 0. 

The difference between the logarithmic price of the European call option 
described by the Ornstein-Uhlenbeck process and the logarithmic price of 
the call option described by the Wiener-Bachelier process is illustrated in 
the Fig. 2. The difference between these logarithmic prices increases with the 




Figure 2: We assume that: k = 0.4, a = 0.03, r = 0, and Cm is the option 
price based on the Merton model, and Cy is the option price based on the 
Vasicek model. 

growth of 9 and fades away in the limit of the expiration date (T = 5) of the 
zero-coupon bond. One of the inappropriate properties of the Merton model 
is the constant drift of the short rate. With the constant positive drift the 
short rate is expected to increase in the future. But this is not realistic. Many 
empirical studies show that, the interest rates exhibit mean reversion, in the 
sense that, if an interest rate is high by historical standards, it typically will 
be falling in the near future (analogically if the current interest rate is low). 
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In paper [2] we considered the prices of options on company's stock sup- 
ported by Wiener-Bachelier model and Ornstein-Uhlenbeck model and we in- 
terpreted the Ornstein-Uhlenbeck process in terms of quantum market games 
theory as a non-unitary thermal tactic. We compared the probability den- 
sity functions of both distributions and showed that the Ornstein-Uhlenbeck 
process would give the true reflection of the market in time scale where the 
approximation by a Wiener-Bachelier process is not valid, but such that mar- 
ket expectations and the asymptotic equilibrium state have not been changed. 
We called this time interval mezzo-scale. The differences between a classical 
look and the pricing which is supported by the quantum model; they are 
visible for very liquid financial instruments. The quantum game theory takes 
advantage the Ornstein-Uhlenbeck process for the financial modelling. Quan- 
tum strategies create unique opportunities for making profits during intervals 
shorter than characteristic thresholds for an effective market (Brownian mo- 
tion). The additional possibilities offered by quantum strategies can lead to 
more successful outcomes than purely classical ones, see [H]. 
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5 Appendix 



In this section we introduce definitions and theorems (without proofs) which 
are used in our paper. They can be find in [3] . [I] . . [TT] . . 

Definition 5.1. 

Let X be an integrable random variable on a probabihty space P), 
and let ^ be a cr-field contained in T. Then the conditional expectation value 
of X, for given is defined to be a random variable E{X\Q) such that: 

• E{X\Q) is ^-measurable, 

• for any A & Q 

I xdp= [ E{x\g)dP. 

J A J A 

Lemma 5.2. 

If Q <Z T is a-field and / : M x — > W is a Borel and bounded function, X 
is a random variable which is independent of Q, and Y >{) is Q -measurable 
and g{y) = E{f{X, y)), then 

E{f{X,Y)\g)=g{Y). 

Theorem 5.3 (Ito lemma). 

(i) Let X = {Xt) t>o be a real-valued process with dynamics 

dXt = atdt + btdWt, 

where a, b are real-valued processes, and W is a one- dimensional stan- 
dard Brownian motion. Then, for any function F : [0, oo) xM 3 [t, x) — » 
F{t, x) e R which is two times continuously differentiable in x and con- 
tinuously differentiable at t, the process defined by Y = {F{t, Xt))t>o is 
an ltd process with dynamics 

/dF dF 1(9^F \ 

,y,= (_(,X,) + a.^«,X,) + -.?^frX,))*+ 

+ 6,^(t,X,)dM.',. 

(a) Let W be a one- dimensional standard Brownian motion, a = {a,P), 
b = (7, 6) and {X, Y) be two-dimensional stochastic processes with dy- 
namics: 

dXt = atdt + PtdWt, 
dYt = itdt + 5tdWt. 
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Then, for any function F : [0, +00) x 9 {t, x, y) — > F{t, x,y) G M 
which is two times continuously differentiable in x and y and contin- 
uously differentiate at t, the process Z — {F{t,Xt,Yt))t>o is an ltd 
process with dynamics 



dZt^ 



dF dF dF 

— {t, Xt, Yt) + at—{t, X„ Yt) + ^t—{t, Xt, Yt) + 



d'^F 
^t^it,Xt,Yt) 



dF dF 
dt + { PtQ^it, Xt, Yt) + 6t-^{t, Xt, Yt) ] dWf 



Definition 5.4. 

A stochastic process X — {Xt)t>o is said to be the geometric Brownian 
motion if it is a solution of the stochastic differential equation 

dXt = tx{t)Xtdt + a{t)XtdWt, 

where yu, a : [0, +00) — > R are deterministic functions of time. The initial 
value for process is assumed to be positive, Xq > 0. 

The function /i is the growing rate (drift) and o" is a volatility rate of the 
process X. 

Theorem 5.5. 

Let X be the geometric Brownian motion with /i = 0, then for any t > the 
following conditions are fulfilled: 

(t) Xt > 0, 

(a) a random variable In ^ has a normal distribution, 

(Hi) if a is a volatility rate of the process X, then 

ft ft 

ft 



Xt = Xsexp [ (-i / a{ufdu+ I a{u)dWu ] , 



Var ^In = J a{uYdu. 
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Lemma 5.6. 

IfY — e^ and X ^ N{m, s"^), then for any constant K > 



£;[max(y - K, 0)] = e"'+2« N 
= E{Y)N 



m 



s 

InK 



+ s 



+ s 



-KN 
KN 



m 



m 



s 

InK 



Theorem 5.7. Let T be the expiration date and K the exercise price of the 
European option, and Ht he the price of the underlying assets at the expiration 
date, then the price of the option at the time t <T is equal to 

Ct = BjE''^ [max(/iT - 0) \Tt\ . (14) 

The basic tool for pricing the financial models is the Feynman-Kac theo- 
rem. We assume that x = {xt)t>o is a diffusion process with dynamics given 
by the stochastic differential equation 

dxt = a{xt, t)dt + (5{xt, t)dW*. (15) 

Moreover, we assume that a process r (r^ = r{xt,t)) describes the short 
interest rate and the risk-neutral probability measure exists. 

Let us consider a security with a single payment of Ht = H{xt, T) and 
the expiration date T. Let Vf = V{xt,t) denotes the price of the security at 
the moment t. It is proved that 



V,^E' 



e-!l<^^^')'^'H{xT,T)\Tt 



(16) 



Theorem 5.8 (Feynman-Kac theorem). 

The function V defined by 



V{x, t) 



E' 



e-St^^(^^''^d'H{xT,T)\J^t]) {x 



satisfies the partial differential equation 
dV dV 1 od'^V 

— {x, t) + aix, t)^ix, t) + -(3{x, t) -Q^ix, t) - r{x, t)V{x, t) = 0, (17) 

together with the terminal condition y(x,T) = H{x,T). The process Vf = 
V{xt,t) describes the price of the considered security. The process {Vt)t>o is 
a diffusion process under P* and W* with the drift 



and the volatility rate 



fi{x, t) = r{x, t)V{x, t) 



dV 

cr{x,t) = —{x,t)P{x,t). 



(18) 



(19) 
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